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' Abstract 

| The effect of magnetic fluctuations on the free energy, the order parameter profile 

and the latent heat at the equilibrium point of the first order phase transition to 
^ superconducting state in thin films of type I superconductors is considered. The 

possibility for an experimental observation of the fluctuation change of the order of 
the superconducting phase transition is briefly discussed. Numerical data for the 
order parameter jump and the latent heat of Al films are presented for needs of 
experimental studies. 

G 

q I 1. Introduction 

The investigation of the fluctuation-induced weakly-first order phase transition in type I 
superconductors in a zero magnetic field known as Halperin-Lubensky-Ma (HLM) effect |T] 
has been recently extended to the case of thin superconducting films It has been 
shown j2] that the HLM effect in quasi-two-dimensional (quasi-2D) films is much stronger 
than in bulk (3D) systems 0] . This result opens an opportunity for an experimental 
observation of the effect in suitably chosen superconducting films. 

In this letter we present new and more precise results for the behaviour of the free energy, 
the order parameter, and the latent heat at the equilibrium phase transition temperature 
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for thin films of type I superconductors. The numerical values of physical quantities of 
experimental interest are calculated for Al films. Our results are compared to those given 
in Refs. 12 El H] • A detailed information about the HLM effect of a fluctuation change 
of the order of the superconducting phase transition in a zero mean magnetic field due to 
persisting magnetic fluctuations and the methods of investigation of this phenomenon are 
published in Refs. [UE1IE]. We shall follow the notations from Ref. [Zj for the parameters 
of the Ginzburg-Landau (GL) free energy of superconductors. 

2. Effective free energy 

The starting point of our consideration is the effective free energy density f(ip) = F(if>)/V 
of a type I superconductor with a volume V which is a function of the mean (uniform) 
superconducting order parameter if) =< if>(x) > of the form (see, e.g., Ref. j2]): 

fW = foW + tfW, (i) 

where 

/o(V0 = a|V>| 2 + ^| 4 , (2) 

and 

1 A 
6f(iP) = ~(D-l)k B T^2\n 

k 

In Eqs. (2) - (3), p(ipo) = Po|^| 2 , where po = (Sue 2 / mc 2 ) , a = a^(T — T c0 ) and b > are 
the usual Landau parameters. They are related to the zero temperature coherence length 
£o = (^ 2 /4ma T c0 ) 1//2 , the zero-temperature critical magnetic field H c0 = a T c0 (4:7r/b) 1 ^ 2 , 
and the initial (unrenormalized) critical temperature T c o f?|. The term 5f(ip) in f(if)) 
describes the effect of the magnetic fluctuations. In Eq. (3), the sum over the wave vector 
k is truncated by the upper cutoff A > k = \k\. For films of thickness L and volume 
V = (L LiL 2 ) we shall assume periodic boundary conditions, and a <^ L < A -1 L a ; 
a = 1,2. For quasi-2D film, i.e. films obeying the condition a L < A" 1 , the sum 
in Eq. (1) contains only terms with zero component ko = 2irno/ Lo(rio = 0) of the wave 
vector k = (ko,ki,k2). This means that the thickness Lq should be smaller than the 
magnetic penetration depth A(T) = (A /|^o| 1,/2 ) which gives the characteristic length of 
the magnetic fluctuations; A(0) = A = (&/po a o^co) is the zero-temperature penetration 
depth, and t = (T — T c0 )/T c0 . We choose the cutoff A = (7r/£ ), so we shall study films 
of size L Q < £o- This is consistent with the general requirement £o < A(T) for the validity 
of the GL free energy for the type I superconductors (see, e.g., Ref. [7j). 

For quasi-2D systems, the continuum limit applied to the sum in Eq. (3) yields a simple 
integral over k = (0, k%, A^). Solving this integral we obtain 

w-^t^M+^-^W]- (4) 

The second term in the r.h.s. of Eq. (4) is nonanalytical and cannot be expanded in 
powers of if). An incomplete Landau expansion of this free energy in positive powers of 
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Wo 
k 2 



(3) 
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\ip 1 2 can be performed, provided 

(polVf/A 2 ) « 1 ■ (5) 

This inequality should be satisfied in the stable (Meissner) phase ip{T) > 0, i.e., at the 
absolute minimum f[i/)(T)] of the function f(ip) below the phase transition temperature. 

Hereafter we shall denote by ip{T) the equilibrium value of if), which describes a stable 
(or metastable) phase and is a solution of the equation of state {df /dip) = 0. The 
condition (5) and the problems included in the reminder of this article can be considered 
in terms of auxiliary parameters corresponding to the theory, in which the 5 f (ip)-pait 
of the free energy is neglected, i.e., corresponding to the free energy /o(/0) given by 
Eq. (2). Within this simplified theory the modulus of the equilibrium zero-temperature 
order parameter is given by \ipo\ = IVK^ = 0)1 = ( a oT c o/by^ 2 . This quantity can be used to 
define the reduced order parameter tp = (\ip\/\ip \) > of the general theory represented 
by the complete free energy f(if>). At equilibrium, i.e., ip = ip{T), the condition (5) 
can be written in the form [ip(T) / XqA] < 1. If we suppose that the difference between 
the values of <p(T) for the simplified and complete theories can be ignored, we shall 
obtain that A(T)A > 1, which is consistent within the superconductivity theory; see 
our discussion below Eq. (3). We have recently investigated Eq. (4) with the help of 
the Landau expansion in Ref. |2|. Here we shall consider the free energy f(ip) given by 
Eqs. (1), (2), and (4) without such an expansion. 

Using the reduced order parameter <p the free energy f(ip) can be represented in the form 

f(p) = § {2t p 2 + p 4 + C(l + t ) [(l + /V) In (1 + /V) - /Vln (w 2 )} } , (6) 
where 



and fi = (£oAA ) 2 . 

3. Results and discussion 

It is known [l a that the HLM effect is stronger in type I superconductors with relatively 
small GL parameter k = (Ao/£o)- This is obvious from Eq. (6). That is why we can 
choose as adequate examples the aluminium (AL, k ~ 0.01, £o = 1.6/im, T c w 1.19K, 
H c {0) w 990e) and tungsten (W, k ~ 0.001, £ = 37/iin, T c w 5 mK, ^(0) w 1.15 Oe). 
One cannot be certain about the best choice of the substance for an experimental check 
of the HLM effect before a comprehensive consideration of the experimental problem. 

Our concrete aim is to establish the size of the effect and we shall do the calculations for 
Al, which is probably one of the best candidates for experiments. The numerical values 
of the parameters C = CL and where A = (hc/2^/2eH c0 ^ ), can be calculated from 
the experimental data for T c , H c (0), and £o of Al, given above. Note, that the available 
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Figure 1: Curves representing the free energy (6) for five values of t : t = —0.001 (see 
□-line), t = -0.00127 (+), t = -0.00137 (o), t = -0.001473 (o), t = -0.0018 (— ). 

experimental data vary depending on the way of preparation of samples and the type of 
measurement. Besides, the experimental values used here are for bulk monocrystals of Al 
and differ within 10 - 20% from those for thin Al films of thickness much smaller than 
the value Lq = 0.1/mi considered below. However, the variation of the experimental data, 
used here does not essentially influence our numerical results. 

The free energy density f{<p) for Al films of thickness L = 0.1/im and several values 
of the parameter to is calculated from Eq. (4). The result is shown in Fig. 1. The 
reduced order parameter corresponding to metastable and stable superconducting states 
is calculated with the help of the equations f((p) = 0, (df /dip) = 0, and the stability 
condition d 2 f/dip 2 > (see Fig. 2). 

Fig. 1 exhibits a well established phase transition of first order. This confirms the preced- 
ing results j2] obtained with the help of Landau expansion; note that in Ref. j2| another 
normalization of the effective free energy has been used, which seems to be less convenient 
than the present one. The positive minima of the free energy describe metastable super- 
conducting states (overheated superconductivity). The respective metastable values of 
the reduced order parameter <p are shown by the dashed curve in Fig. 2. The metastabil- 
ity states are closed in the temperature interval from T* = 0.9990T c0 to T eq « 0.9985T c0 
- the equilibrium temperature of the first order phase transition. The curve marked by 
squares (□) in Fig. 1 corresponds to T = T*, i.e., this is the curve, at which the minimum 
of the free energy for a nonzero value of ip appears for the first time when the tempera- 
ture is lowered from the side of the normal phase. The curve drawn by " diamonds" (o) 
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Figure 2: Order parameter profile near T c0 . The vertical line at t — —0.001473 indicates 
the equilibrium jump of the order parameter. 



corresponds to the equilibrium transition temperature T eq , at which the minimum of the 
free energy for nonzero value of the order parameter (p is equal to zero. There the energies 
of the superconducting and the normal phases are equal. The superconducting phase is 
stable for all temperatures below T eq . The stable superconducting states are shown by 
the solid curve in Fig. 2 for several values of t below t 0eq = t (T eq ) m —0.0015. 

The equilibrium entropy jump 5s = —df[T,ip(T)]/dT at T eq per unit volume (the total 
entropy jump is 5S = V5s) is obtained in the form 



5s 



Hi 



AtxT, 



c0 



c 



eq 



[(i 



Ig) 111 (1 + Wlq) - Wlq^ {wlq)] 



(8) 



For the Al film of thickness Lq = 0.1/j.m, the C— term in the curly brackets of Eq. (8) 
is of order 10~ 3 ipl q and can be neglected. Taking the value of the order parameter jump 
(p eq ~ 0.033 from Fig. 2, we obtain 



5s = 



tt2 

co -v? 2 



c0 



eq 



-0.7 



erg 
K.cm 3 



(9) 



In order to show the significance of this result we compare it with the specific heat capacity 
jump at the second order phase transition point, T c0 , described by the free energy fo(p)'- 



SC(T c0 ) = 



H c0 

4nT r 



660 



c0 



erg 
K.cm 3 



(10) 
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The ratio (Ss/5C) ~ 0.001 is 10 3 times bigger than the respective quantity for bulk Al PQ. 
Comparing with the results [3] for bulk Al one can see that the metastability domain 
(T* — T eq ) for the film under consideration is also 10 3 times larger than for bulk Al. The 
calculated value ((p eq ~ 0.0032) of the order parameter jump in bulk Al is about 0.1 of 
that shown in Fig. 2. While the overheated and stable superconducting states in quasi-2D 
Al films appear only for to < 0, i.e. below T co , in 3D Al these states occur also slightly 
above T c0 . 

The present results are consistent with those obtained by Landau expansion [2] of the 
free energy (6). However, there is a difference in the numerical values of some parameters 
considered in the present investigation and those in Ref. Obviously, the Landau 
expansion has some limitations when applied to free energies of the type (6). 

4. Conclusion 

The results show that Al films of thickness below 0.1/im can be used for an experimental 
test of the HLM effect. Calculations for thinner Al films or for W, where the GL parameter 
is 10 times smaller, will give much stronger effect and, hence, better opportunity for an 
experimental test. The thickness of the film may be lowered up to the nm scale, below 
which the superconductivity is destroyed. 
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